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Abstract. It is generally assumed that as a result of multiple scattering, the source 
CN| , distribution measured in HBT interferometry corresponds to a chaotic source at 

freeze-out. This assumption is subject to question as effects of multiple scattering 
^ ' in HBT measurements must be investigated within a quantum-mechanical framework. 

XjT^ , Applying the Glauber multiple scattering theory at high energies and the optical model 

at lower energies, we find that multiple scattering leads to an effective HBT density 
. distribution that depends on the initial chaotic source distribution with an absorption. 

o ■ 
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1. Introduction 



p • Recent experimental measurements of HBT correlations in relativistic heavy-ion 
p • collisions show only relatively small changes of the extracted longitudinal and transverse 
^ ' radii as a function of collision energies, and the ratio of i?out/-Rside ~ 1 [!]■ The 
difficulties of explaining these HBT measurements with theoretical models, known as the 
' "HBT puzzles", have been discussed by many authors [2, 3]. In these comparisons with 

;h ' theoretical models, it is generally assumed that as a result of multiple scattering, the 
source distribution measured in an HBT measurement corresponds to a chaotic source 
at freeze-out, in which a detected hadron suffers its last hadron-hadron scattering. 

In a recent quantum-mechanical treatment of the multiple scattering process using 
the Glauber theory at high energies and the optical model at lower energies, it was found 
that the HBT interferometry does not measure the freeze-out source distribution and 
the effective HBT density distribution depends on the initial chaotic source distribution 
with an absorption [4]. Effects of collective flows have also been investigated [4]. 

What is the physical basis for these new insights in HBT measurements? While 
the detailed arguments leading to the above results have been presented previously 
in Ref. [4], it is instructive to review here the most important features of the multiple 
scattering process and the HBT interferometry that give rise to the above unconventional 
viewpoints. We shall first examine the origin of the HBT interferometry in Section II and 
study next how the multiple scattering process will affect the HBT correlation function 
in Section HI, using the quantum- mechanical Glauber theory of multiple scattering [5]. 

At the recent Quark Matter 2004 Conference, Kapusta and Li [6] reported new 
findings which support qualitatively some of the earlier results of Ref. [4]. 
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2. HBT Interferometry for Chaotic Sources 

We examine briefly the origin of the HBT correlation for identical bosons (pions), as 

reviewed in Chapter 17 of Rcf. [7]. The HBT correlation is represented by the probability 
P{ki,k2) of detecting identical bosons with 4-momenta ki and k2 in coincidence, or 
alternatively, by the correlation function R{ki,k2) = [P{ki,k2)/P{ki)P{k2)] -1 where 
P{k) is the single-particle momentum distribution. 
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Figure 1. The shaded region represents the initial source which will 
expand to the freeze-out configuration before the bosons reach the detectors. 
P{ki, k2) is given by squaring the total amplitude from all source points (xi, a;2), 

P{k,,k2) = ^1 E >l(A;iXi)e^'^°('=i-i)>l(A;2X2)e^'^°(*^^-^)Vi2(2:i,X2)|', (1) 

where A{kiXi) and (f)Q{kiXi) are the amplitude and the phase for the production of boson 
ki at the space-time point Xi (Fig. 1). ipi2{xi,X2) is the amplitude for the propagation of 
a pair of bosons {ki,k2) from the source points {xi, X2) to the detection points (x^i, 0:^2), 



V2 



^iki-{xdi-xi)+ik2-(xd2-X2) _j_ ^iki-{xdi-X2)+ik2-(xd2-xi) I. 



where the two terms represent amplitudes for two different histories of traveling from the 
source points to the detection points (the solid and the dashed trajectories (histories) 
of Fig. 1). If the source is coherent, there is no HBT correlation and R{ki, ^2) = 0. 

If the source is chaotic with random and fluctuating phases (f)o{kx), the absolute 
square of the sum in Eq. (1) becomes the sum of absolute squares, 

P{k,,k2) = E A\k,X,)A''{k2X2Ml2{Xi,X2)\^. (3) 



Xl,X2 



The HBT correlation is then present with a non- vanishing R{ki,k2). The correlation 
function -R(/ci, /C2) depends on the cross term of | ■0i2 (3:^1, 2:2) P of Eq. (2), which contains 
the phase difference between the two histories, 

ki ■ {Xdl - Xi) + k2 ■ {Xd2 - X2)- {ki ■ {Xdl - X2) + k2 ■ {Xd2 - Xi)} 

= -{ki - k2) ■ xx - {k2 - ki) ■ X2. (4) 
Representing Y^xi by / dxip{xi) with the source spatial density p(xi), Eqs. (2)- (4) give 

R{kx,k2) -\ Jdx e-'^^'-^^>^ Pes{x;kuk2) \\ (5) 

p{x)A{ki,x)A{k2,x) 
Peff (x; /ci , /C2) = , ■ (6) 



IPiki)Pik2) 

The measurement of R{ki,k2) therefore provides information on the Fourier transform 
of the effective source distribution Pes{x; ki, /C2) of a chaotic source. 
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3. Multiple Scattering and HBT 

In high-energy heavy-ion coUisions, particles such as pions are produced and they will 
undergo multiple scattering until the system reaches the freeze-out configuration (Fig. 
1). According to the Glauber theory of multiple scattering [5], the probability amplitude 
for a pion with momentum ki to go from the source point Xi to the freeze-out point x/i 
and then to the detection point Xdi is 

Xdi) = exp{iki ■ {xdi - xi) + i(t)s{xfi - xi)}. (7) 

At high particle energies, the phase ^^(x/i — xi) is the sum of two-body scattering phase 
shifts Xj of the particles with which the particle scatters along its trajectory, 

Af(a:yi||-a;i||) 

,{xfi -xi)= xA^ii. - Xj±), (8) 
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where Xi\\ and X\± are coordinates longitudinal and transverse to ki. At lower energies, 
'Psixfi — xi) is given in terms of the two-body optical potential Vj{xi — Xj) and the 
particle velocity v, N{xnu-xiii) 

4>s{xfi -xi) = - - ^ Vj{xi - Xj)dxi\\. (9) 



The probability amplitude (7) leads to the proper classical transport description of the 
particle in a medium [8]. It modifies the wave function ■0i2 (3^1)^2) of Eq. (2) to 

1p^2{Xl X2) = — — { Q^^l-{^dl-Xl)+i(j>s(x fl-X^)+ik2-{xd2-X2)+i<t>s{x f2-X2) 

_|_ ^ik\_-(xdi-X2)+i(t>s(,x'j,2-X2)-\-ik2-{xd2-xi)+i<t>a{x'j^-xi) I ^-j^g^ 

where and {x'^ix'^^) are the two sets of freeze-out coordinates for the two 

histories (Fig. 1). In an HBT measurement of a chaotic source, the difference of the 
phases in the cross term of |'0i2(a^i, a;2)P in Eq. (4) is now modified to 

^1 • {Xdl - Xi) (t)s{Xfi - Xi) k2 ■ {Xd2 - X2) + (f>s{Xf2 " X2) 
-{h ■ {Xdl - X2) + (t>*s{x'f2 - X2) + k2 ■ {Xd2 - Xi) + (f)*{x'fi - Xi)} 

= -{ki - k2) ■ xi - {k2 - ki) ■X2 + A, (11) 

where we introduce A to represent the effects of multiple scattering. For the 
measurement of i?out, — ^2 is along ki, we have Xfi — x'f^ and Xf2 — x'^21 

(f)s{xfi- Xi) ^ (psix'fi- xi), and (ps{xf2- X2) = (psix'j:^- X2). (12) 

For the measurement of -Rside and -Riong, fci — ^2 = q is perpendicular to ki, and 

Tie 4>s{xfi - xi) - Tie (psix'fi - Xi) = {xfi - x'j-^) ■ VTZe 0s(a;/i - Xi) 

oc q ■ V ±7le (f)s{xfi — xi) ~ 0, (13) 

where the last approximate equality arises as the vector sum in V±Tle (pg with random 
transverse vector contributions from many independent scatterings is approximately 
zero. In both cases, the real parts of the phase differences cancel approximately, and 
only the imaginargy parts remain, 

A ~ 2i Im 4>s{xfi — x\) -\- 2i Tm (f)s{xf2 — X2). (14) 
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From Eqs. (11) and (14), R{ki, /c2) is again given by Eq. (5) but the HBT effective density 
distribution PeG{x] ki, k2) of Eq. (6) is modified, in tfie case of multiple scattering, to [4] 



The effective distribution depends on the initial source distribution with an absorption. 
4. Conclusions and Discussions 

Comparing the results of the last two sections, we see that the new insights concerning 
the multiple scattering process and HBT measurements arise from the following 
considerations: (1) In a quantum mechanical description, the multiple scattering process 
gives rise to the accumulation of phases along the trajectories of the detected particles. 
(2) For a chaotic source, the HBT correlation arises from the difference in the phases 
accumulated in two different sets of particle trajectories. (3) For these two sets of 
trajectories, the real parts of the accumulated phases due to the multiple scattering 
process approximately cancel each other and only the imaginary absorptive parts remain. 

Based on these considerations, we find that the multiple scattering process leads to 
an effective density distribTition that depends on the initial chaotic source distribution 
with an absorption, Eq. (15). The effects of the longitudinal momentum loss can be 
further included in the future by using the extension of the Glauber theory formulated 
by Blankenbecler and Drell [9]. While the absorption and the longitudinal momentum 
loss will modify the transmission of the initial chaotic source distribution, the effective 
source is closer to the initial chaotic source configuration and the nuclear geometrical 
overlap than the freeze-out configuration. As a consequence, the present new insights 
may pave the way for a better understanding of the HBT puzzles. 
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